The impact of sequencing depth on the inferred taxonomic composition and AMR gene content of metagenomic samples

Shotgun metagenomics is increasingly used to characterise microbial communities, particularly for the investigation of antimicrobial resistance (AMR) in different animal and environmental contexts. There are many different approaches for inferring the taxonomic composition and AMR gene content of complex community samples from shotgun metagenomic data, but there has been little work establishing the optimum sequencing depth, data processing and analysis methods for these samples. In this study we used shotgun metagenomics and sequencing of cultured isolates from the same samples to address these issues. We sampled three potential environmental AMR gene reservoirs (pig caeca, river sediment, effluent) and sequenced samples with shotgun metagenomics at high depth (~ 200 million reads per sample). Alongside this, we cultured single-colony isolates of Enterobacteriaceae from the same samples and used hybrid sequencing (short- and long-reads) to create high- quality assemblies for comparison to the metagenomic data. To automate data processing, we developed an open- source software pipeline, ‘ResPipe’

Genetic, epigenetic and microbiome characterisation of an earthworm species (Octolasion lacteum) along a radiation exposure gradient at Chernobyl

The effects of exposure to different levels of ionising radiation were assessed on the genetic, epigenetic and microbiome characteristics of the “hologenome” of earthworms collected at sites within the Chernobyl exclusion zone (CEZ). The earthworms Aporrectodea caliginosa (Savigny, 1826) and Octolasion lacteum (Örley, 1881) were the two species that were most frequently found at visited sites, however, only O. lacteum was present at sufficient number across different exposure levels to enable comparative hologenome analysis. The identification of morphotype O. lacteum as a probable single clade was established using a combination of mitochondrial (cytochrome oxidase I) and nuclear genome (Amplified Fragment Length Polymorphism (AFLP) using MspI loci). No clear site associated differences in population genetic structure was found between populations using the AFLP marker loci. Further, no relationship between ionising radiation exposure levels and the percentage of methylated loci or pattern of distribution of DNA methylation marks was found. Microbiome structure was clearly site dependent, with gut microbiome community structure and diversity being systematically associated with calculated site-specific earthworm dose rates. There was, however, also co-correlation between earthworm dose rates and other soil properties, notably soil pH; a property known to affect soil bacterial community structure. Such co-correlation means that it is not possible to attribute microbiome changes unequivocally to radionuclide exposure. A better understanding of the relationship between radionuclide exposure soil properties and their interactions on bacterial microbiome community response is, therefore, needed to establish whether these the observed microbiome changes are attributed directly to radiation exposure, other soil properties or to an interaction between multiple variables at sites within the CEZ

Citizen science monitoring reveals links between honeybee health, pesticide exposure and seasonal availability of floral resources

We use a national citizen science monitoring scheme to quantify how agricultural intensifcation afects honeybee diet breadth (number of plant species). To do this we used DNA metabarcoding to identify the plants present in 527 honey samples collected in 2019 across Great Britain. The species richness of forage plants was negatively correlated with arable cropping area, although this was only found early in the year when the abundance of fowering plants was more limited. Within intensively farmed areas, honeybee diets were dominated by Brassica crops (including oilseed rape). We demonstrate how the structure and complexity of honeybee foraging relationships with plants is negatively afected by the area of arable crops surrounding hives. Using information collected from the beekeepers on the incidence of an economically damaging bee disease (Deformed Wing Virus) we found that the occurrence of this disease increased where bees foraged in agricultural land where there was a high use of foliar insecticides. Understanding impacts of land use on resource availability is fundamental to assessing long-term viability of pollinator populations. These fndings highlight the importance of supporting temporally timed resources as mitigation strategies to support wider pollinator population viability

Genomic network analysis of environmental and livestock F-type plasmid populations

F-type plasmids are diverse and of great clinical significance, often carrying genes conferring antimicrobial resistance (AMR) such as extended-spectrum β-lactamases, particularly in Enterobacterales. Organising this plasmid diversity is challenging, and current knowledge is largely based on plasmids from clinical settings. Here, we present a network community analysis of a large survey of F-type plasmids from environmental (influent, effluent, and upstream/downstream waterways surrounding wastewater treatment works) and livestock settings. We use a tractable and scalable methodology to examine the relationship between plasmid metadata and network communities. This reveals how niche (sampling compartment and host genera) partition and shape plasmid diversity. We also perform pangenome-style analyses on network communities. We show that such communities define unique combinations of core genes, with limited overlap. Building plasmid phylogenies based on alignments of these core genes, we demonstrate that plasmid accessory function is closely linked to core gene content. Taken together, our results suggest that stable F-type plasmid backbone structures can persist in environmental settings while allowing dramatic variation in accessory gene content that may be linked to niche adaptation. The association of F-type plasmids with AMR likely reflects their suitability for rapid niche adaptation

Niche and local geography shape the pangenome of wastewater- and livestock-associated Enterobacteriaceae

Escherichia coli and other Enterobacteriaceae are diverse species with “open” pangenomes, where genes move intra- and interspecies via horizontal gene transfer. However, most analyses focus on clinical isolates. The pangenome dynamics of natural populations remain understudied, despite their suggested role as reservoirs for antimicrobial resistance (AMR) genes. Here, we analyze near-complete genomes for 827 Enterobacteriaceae (553 Escherichia and 274 non-Escherichia spp.) with 2292 circularized plasmids in total, collected from 19 locations (livestock farms and wastewater treatment works in the United Kingdom) within a 30-km radius at three time points over a year. We find different dynamics for chromosomal and plasmid-borne genes. Plasmids have a higher burden of AMR genes and insertion sequences, and AMR-gene-carrying plasmids show evidence of being under stronger selective pressure. Environmental niche and local geography both play a role in shaping plasmid dynamics. Our results highlight the importance of local strategies for controlling the spread of AMR